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Summary  of  Research  Progress 


Objective 

Catalysis  of  heat-releasing  H+OH  recombination  in  the  nozzle  section  of  a  SCRAM-jet 
potentially  increases  the  fuel  efficiency  and  therefore  the  thrust  of  a  hydrogen-fueled,  air-breathing, 
hypersonic  aircraft.  Phosphine  (PH3)  combustion  products  substantially  increase  the  rate  of 
disappearance  of  excess  OH  when  present  in  a  mixture  of  partially  photolyzed  water  vapor.  The 
objectives  of  this  research  program  were  to  1)  determine  the  identity  and  equilibrium  densities  of 
phosphorus  combustion  products  present  in  the  gas  phase  when  phosphine  is  burned  with  oxygen, 
2)  construct  a  reaction  model  including  the  most  abundant  phosphorus  species,  3)  determine  the 
rate  constants  for  reactions  of  phosphorus  combustion  products  in  a  high  temperature  gas,  and 
4)  identify  the  mechanism  by  which  phosphorus  combustion  products  increase  the  disappearance 
rate  of  excess  OH  density  in  a  water  photolysis  experiment.  To  assess  the  potential  effectiveness  of 
H+OH  ^  H2O  catalysis,  the  reaction  mechanism  developed  in  this  study  was  used  in  a  computer 
model  of  SCRAM-jet  reaction  kinetics  to  predict  the  gain  in  heat  release  in  the  nozzle  from 
phosphorus  combustion  products  and  to  calculate  the  increase  in  specific  impulse  resulting  from 
the  heat  release  gain. 

Results 

Photolysis  Experiments 

Hydrogen  reacts  with  oxygen  at  high  temperature  to  produce  water  vapor  and  a 
considerable  amount  of  heat.  In  our  combustion  reactor,  gas  temperatures  of  1500-2500  K  are 
achieved  at  pressures  of  300-700  Ton*.  At  these  temperatures  the  equilibrated  combustion  products 
consist  of  small  but  significant  amounts  of  radicals  such  as  OH  and  H.  To  probe  the  reaction 
chemistry  of  the  combustion  products,  the  equilibrated  system  is  perturbed  by  laser  photolysis  of 
H2O,  which  substantially  increases  the  pool  of  OH  and  H  radicals. 

In  one  set  of  experiments  [1],  the  return  of  the  system  to  equilibrium  is  experimentally 
investigated  by  recording,  as  a  function  of  time,  the  OH  absorption  of  a  frequency-doubled  cw-dye 
laser  probe  beam  producing  312  nm  radiation.  The  OH  absorption  signal  is  recorded  on  a  digital 
oscilloscope  and  transferred  as  a  binary  file  to  a  computer  for  subsequent  analysis. 
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Figure  1  shows  the  dramatic  effect  on  the  rate  of  recovery  of  the  OH  signal  when  PH3  is 
added  to  the  hydrogen  fuel  and  burned  with  oxygen.  Photolysis  occurs  at  0.1  ms  on  the  graph  and 
results  in  a  very  fast  increase  in  OH  absorbance  (the  value  plotted  in  Fig.  1  is  the  change  in  OH 
absorbance  from  the  equilibrium  level).  When  phosphine  is  not  present  in  the  hydrogen  fuel 
stream,  the  OH  signal  relaxes  back  towards  its  equilibrium  value  at  a  slow  rate  represented  by  the 
empty  squares  plotted  in  Fig.  1.  The  rate  at  which  OH  absorbance  returns  to  its  equilibrium  value 
is  increased  dramatically  by  the  addition  of  0.1%  by  volume  phosphine  in  the  hydrogen  flow 
stream  (filled  circles).  The  addition  of  phosphine  in  amounts  as  small  as  0.03%  by  volume  (not 
shown  in  this  figure)  results  in  a  significant  increase  in  OH  disappearance  rate. 


Time  (microseconds) 

Fig.  1  OH  density  change  following  H2O  photolysis.  Test  gas  mixture  is  the  combustion  products 
of  H2(empty  squares)  and  0.1%  PH3/H2  (filled  circles).  Rich  fuel  mixture. 


The  more  rapid  disappearance  of  excess  OH  when  PH3  combustion  products  are  present  in 
the  photolyzed  test  gas  is  not  proof  of  an  increased  rate  of  H+OH  recombination.  One  may  argue 
that  a  fast  bimolecular  reaction  step  consumes  OH  but  does  not  generate  H2O  until  much  later.  To 
more  completely  test  the  hypothesis  that  PH3  combustion  products  catalyze  the  overall  rate  of 
H+OH  recombination,  the  product  side  of  this  reaction  was  measured  [2].  An  H2O  absorbance 
signal  at  192  nm  resulted  from  a  rise  in  temperature  caused  by  heat  released  when  the  photolysis 
products,  H  and  OH,  recombine. 

Figure  2  shows  a  typical  H2O  absorbance  signal  for  a  fuel-rich  test  gas  mixture.  The 
absorbance  change  shown  in  the  figure  was  calculated  by  taking  the  decadic  logarithm  of  the 
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average  signal  (time  series  of  photodetector  output)  recorded  prior  to  the  photolysis  pulse  divided 
by  the  signal  itself.  The  signal  slowly  increases  to  a  positive  limiting  value  of  about  0.003  after 
photolysis.  The  shape  of  the  signal  suggests  the  appearance  of  a  product  developing  over  the 
0.9  ms  experimental  time  period.  For  this  experiment,  photolysis  of  H2O  by  the  excimer  laser 
creates  excess  OH  and  H  radical.  When  the  radicals  eventually  recombine  the  laser  energy  used  in 
photolyzing  H2O  appears  as  thermal  energy  which  raises  the  temperature  of  the  gas.  The  signal  at 
192  nm  is  in  a  spectral  region  where  the  absorbance  of  H2O  increases  with  temperature.  Therefore, 
the  transient  absorbance  signal  at  192  nm  has  two  components  that  are  superposed,  an  initial 
decrease  in  the  absorbance  signal  due  to  a  decrease  in  H2O  density  and  a  gradual  increase  in  the 
absorbance  signal  due  to  a  change  in  the  absorption  coefficient  caused  by  an  increase  in 
temperature.  Calculations  show  that  the  signal  expected  at  192  nm  is  predominately  determined  by 
the  temperature-driven  increase  in  the  absorption  coefficient  of  H2O  when  recombination  reactions 
release  the  photolysis  energy  stored  in  the  population  of  excess  OH  and  H  radicals. 


Fig.  2  H2O  absorbance  change  following  H2O  photolysis.  Test  gas  mixture  is  the  combustion 
products  of  H2(empty  squares)  and  0.1%  PH3/H2(filled  circles).  Rich  fuel  mixture. 

As  Fig.  1  shows,  phosphine  combustion  products,  when  present  at  relatively  low  density, 
substantially  increase  the  rate  at  which  the  OH  density  returns  to  its  equilibrium  value  following 
H2O  photolysis.  If  the  more  rapid  relaxation  rate  for  the  OH  signal  is  the  result  of  an  increased  net 
H+OH  recombination  rate,  then  the  H2O  signal  should  rise  more  quickly  to  its  limiting  value  when 
PH3  combustion  products  are  present  in  the  test  gas  mixture.  Figure  2  compares  the  H2O 
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absorbance  signals  recorded  with  and  without  0.1%  PH3  present  in  the  H2  fuel  stream  for  a  fuel 
rich  combustion  mixture.  A  faster  rise  in  the  H2O  signal  at  earlier  times  is  observed  when  PH3  is 
added  to  the  fuel  indicating  a  faster  release  of  heat  caused  by  more  rapid  H+OH  recombination. 

The  observed  H2O  signals  are  consistent  with  the  OH  measurements  and  model 
calculations. 

Reaction  Model 

A  kinetic  model  was  constructed  to  analyze  the  experimental  OH  and  H2O  absorbance 
signals  and  the  effect  of  phosphine  combustion  products  on  these  signals.  A  list  of  those 
phosphorus  species  (PO,  PO2,  PO3,  HOPO,  HOPO2,  HPO,  P2O3,  H2POH,  HPOH,  P2O2,  P2O, 
PH3,  PH2,  PH,  P4,  P2,  P)  that  potentially  participate  in  the  reaction  chemistry  was  assembled,  and 
the  thermodynamic  properties  of  each  species  was  either  taken  from  literature  reports  or  estimated 
from  quantum  chemical  calculations  reported  in  the  literature  [1]. 

In  the  original  reaction  model  proposed  [1],  the  equilibrium  densities  of  the  17  phosphorus 
species  were  calculated  under  conditions  of  the  photolysis  experiment.  Species  that  did  not  have  a 
large  enough  density  to  cause  an  increase  in  the  H+OH  recombination  rate  even  at  collision- 
controlled  rates  were  ehminated  from  further  consideration.  A  reaction  model  was  constructed  from 
reactions  of  the  remaining  species. 

In  subsequent  work  [3],  a  computer  program  was  written  to  list  all  possible  bimolecular 
and  recombination  reactions  occurring  among  the  17  phosphorus  species.  Another  program 
estimated  the  rate  constants  for  the  bimolecular  reactions;  rate  constants  for  the  recombination 
reactions  were  estimated  using  a  simplified  RRKM  calculation  [1].  A  reaction  model  reduction 
analysis  was  used  to  select  those  reactions  that  contribute  to  the  overall  rate  of  H+OH 
recombination  in  a  gas  that  has  been  initially  removed  from  chemical  equilibrium  by  H2O 
photolysis  or  by  a  temperature  jump  perturbation.  A  thermodynamic  free  energy  function  was  used 
to  quantify  the  rate  of  system  relaxation  back  to  equiUbrium  for  a  series  of  36  reaction  conditions 
covering  a  temperature  range  of  1500  to  3000  K,  a  gas  density  range  of  5  x  10'^  to  5  x  10'^ 
moles/cm^  and  a  fuel  equivalence  ratio  of  0.8  to  1.2.  The  reactions  that  show  the  greatest 
participation  in  equilibrium  relaxation  kinetics  by  this  analysis  are  listed  below  in  Table  1.  Criterion 
for  inclusion  of  a  reaction  in  the  table  is  that  the  reaction  have  an  analysis  factor  (defined  in 
reference  3)  larger  than  0.0001  under  any  set  of  gas  conditions  for  either  temperature  or  photolysis 
perturbation. 
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Table  1.  Selected  Reactions  and  Associated  Analysis  Factors 


Reaction 

Analysis  Factor 

Temp 

Photo 

OH  +  PO  =  H  +  P02 

0.0017 

0.0009 

P02  +  P  =  PO  +  PO 

0.0015 

0.0013 

0  +  P03  =  02  +  P02 

0.0013 

0.0009 

H20  +  P03  =  OH  +  H0P02 

0.0009 

0.0027 

0  +  P02  =  P03 

0.0006 

0.0083 

OH  +  PO  =  HOPO 

0.0005 

0.0001 

OH  +  P  =  H  +  PO 

0.0004 

0.0004 

H  +  P203  =  PO  +  HOPO 

0.0003 

0.0097 

OH  +  HOPO  =  H  +  H0P02 

0.0003 

0.0008 

02  +  P  =  0  +  PO 

0.0003 

0.0003 

0  +  HPO  =  OH  +  PO 

0.0002 

0.0002 

H  +  P03  =  OH  +  P02 

0.0002 

0.0013 

H0P02  +  P  =  PO  +  HOPO 

0.0002 

0.0000 

H  +  HPO  =  HPOH 

0.0001 

0.0002 

PO  +  HPOH  =  HOPO  +  PH 

0.0001 

0.0000 

0  +  P203  =  P02  +  P02 

0.0001 

0.0003 

H  +  P03  =  0  +  HOPO 

0.0001 

0.0010 

0  +  P203  =  PO  +  P03 

0.0000 

0.0027 

0  +  P2  =  PO  +  P 

0.0000 

0.0001 

Reaction 

Analysis  Factor 

Temp 

Photo 

H  +  P02=HOPO 

0.7958 

0.6878 

H  +  HOPO  =  H2  +  P02 

0.5104 

0.5136 

OH  +  P02  =  HOP02 

0.4139 

0.2911 

PO  +  P02  =  P203 

0.3517 

0.0009 

0  +  HOP02  =  02  +  HOPO 

0.1437 

0.0710 

0  +  HOPO  =  OH  +  P02 

0.0830 

0.0433 

OH  +  HOPO  =  H20  +  P02 

0.0821 

0.0054 

PO  +  HOP02  =  P02  +  HOPO 

0.0791 

0.0002 

H  +  HOPO  =  H20  +  PO 

0.0762 

0.0106 

H  +  HOP02  =  H20  +  P02 

0.0329 

0.0170 

0  +  HOPO  =  HOP02 

0.0233 

0.0087 

0  +  PO  =  P02 

0.0083 

0.0038 

H02  +  P02  =  OH  +  P03 

0.0054 

0.0005 

H02  +  P02  =  0  +  HOP02 

0.0053 

0.0004 

02  +  PO  =  0  +  P02 

0.0051 

0.0045 

P2  +  P2  =  P4 

0.0038 

0.0000 

H02  +  P02  =  02  +  HOPO 

0.0038 

0.0001 

H  +  PO  =  HPO 

0.0027 

0.0028 

H  +  HPO  =  H2  +  PO 

0.0021 

0.0021 

OH  +  P203  =  P02  +  HOPO 

0.0018 

0.0000 

Most  of  the  reactions  originally  selected  for  inclusion  in  a  reaction  model  (in  reference  1) 
are  present  in  Table  1.  Conclusions  reached  concerning  catalysis  reaction  schemes  that  are  of 
importance  for  H+OH  recombination  are  unchanged  by  the  more  rigorous  analysis.  The  most 
likely  reaction  sequence  catalyzing  both  H+H  and  H+OH  recombination  is 

H  +  PO2  ^  HOPO 
H  +  HOPO  -+  H2  +  PO2 
OH  +  H2  ->  H2O  +  H 
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Temperature  Dependence  of  Rate  Constants  [5] 

The  relaxation  rate  of  OH  density  back  to  equilibrium  following  H2O  photolysis  was 
recorded  under  fuel  rich  conditions  at  gas  temperatures  in  the  range  of  1500-2500  K.  A  reaction 
model  consisting  of  a  standard  set  of  13  H2/O2  reactions  and  seven  of  the  most  important  (as 
determined  by  the  reaction  model  reduction  analysis)  phosphorus  reactions  was  used  to  fit  the  OH 
decay  signals. 

The  laser  photolysis  experiments  consisted  of  two  sets  of  data;  those  in  which  about  0.1% 
phosphine  in  hydrogen  was  burned  with  oxygen  and  those  in  which  phosphine  was  not  included. 

The  pure  H2  data  set  was  analyzed  by  allowing  the  rate  constants  for  all  13  reactions  to 
vary.  After  the  optimal  set  of  rate  constants  was  calculated  (those  that  minimized  the  deviation  of 
computed  from  experimental  OH  absorbance  signal),  each  rate  constant  was  allowed  to  vary 
separately  and  a  range  of  values  tor  each  rate  constant  was  calculated  over  which  the  deviation  did 
not  change  by  more  than  10%.  Most  of  the  rate  constants  had  a  very  wide  range  of  values 
indicating  the  relative  insensitivity  of  the  computed  OH  signal  to  the  value  of  the  rate  constant. 
However,  the  deviation  of  computed  from  experimental  OH  absorbance  signal  showed  the 
expected  relatively  high  sensitivity  to  the  value  of  the  rate  constant  for  the  H  +  OH  +  M  H2O  + 
M  reaction.  The  rate  constant  for  this  recombination  reaction  displayed  a  temperature  dependence 
that  differed  substantially  from  that  of  the  accepted  literature  value.  Figure  3  shows  a  graph  of 
Logio[k/ko]  versus  temperature  for  the  H+OH  recombination  reaction,  where  k  is  rate  constant 
derived  from  the  analysis  and  ko  =  8.32  x  lO^^T'^. 

The  PH3/H2  data  was  analyzed  in  a  manner  similar  to  that  of  the  pure  H2  data.  However, 
for  this  data  set  only  the  seven  rate  constants  associated  with  the  seven  phosphorus-containing 
reactions  were  used  as  parameters  in  the  fit  (the  H2/O2  rate  constants  were  fixed  at  their  standard 
values).  The  result  of  this  analysis  showed  an  insensitivity  of  the  computed  OH  signal  to  all  but  the 
H  +  PO2  HOPO  and  H  +  HOPO  ^  H2  +  PO2  reactions.  While  both  these  reactions  showed  a 
temperature  dependence  that  deviated  at  some  temperatures  from  the  baseline  values,  a  trend  was 
not  evident  in  the  data.  In  Fig.  4  below,  a  graph  of  Logio[k/ko]  versus  temperature  for  the  H  + 
PO2  HOPO  reaction  is  shown. 
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Fig.  4  Logjolk^]  versus  temperature  where  k  is  the  rate  constant  calculated  for  the  H  +  PO2  — > 
HOPO  reaction  and  ko  is  the  baseline  value  of  this  rate  constant. 
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Nozzle  Performance  Calculations  [4] 

It  is  generally  believed  that  hypersonic  aircraft  will  suffer  from  diminished  fuel  efficiency 
as  a  consequence  of  slow  radical  recombination  in  the  aircraft  nozzle.  A  faster  recombination  rate 
may  allow  increased  heat  release  in  the  nozzle  and  increased  thrust  when  the  thermal  energy  is 
converted  to  directed  gas  flow  by  nozzle  expansion.  Whether  a  homogeneous  catalyst,  such  as  a 
mixture  of  phosphine  combustion  products,  can  increase  the  overall  system  performance  of  a 
hypersonic  vehicle  will  depend  on  the  trade-off  between  the  extra  weight  of  the  catalyst  and  the 
benefit  realized  by  the  additional  energy  released  in  the  nozzle. 

A  series  of  simple  one-dimensional  nozzle  expansion  calculations  in  which  the  heat  released 
per  unit  of  fuel  mass  was  compared  with  and  without  phosphine  addition  to  the  fuel  were 
performed.  The  calculations  assumed  adiabatic  reaction  of  hydrogen  fuel  with  air.  A  chemically 
equilibrated  gas  of  combustion  products  entered  the  nozzle  section  of  the  SCRAM-jet.  The  nozzle 
was  assumed  to  be  rectangular  in  cross  section  and  expand  in  one  dimension  only.  The  heat 
released  per  unit  fuel  mass  by  the  flowing  gas  in  the  nozzle  was  calculated  as  a  function  of  nozzle 
length.  Comparisons  between  pure  H2  fuel  and  H2  containing  PH3  revealed  that  the  latter  released 
additional  heat  throughout  the  length  of  the  nozzle  for  fuel-rich  combustion  conditions,  as  shown 
in  Fig.  5. 


Nozzle  Distance  (m) 

Fig.  5  Fractional  increase  in  heat  release,  (j),  versus  nozzle  distance.  Fuel  mixture  has  an 
equivalence  ratio  of  nominally  1.2,  and  the  gas  velocity  at  the  nozzle  entrance  is  5000  m/s. 
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The  difference  in  specific  impulse  of  the  pure  H2  and  H2/PH3  fuels  was  also  calculated.  In 
Fig.  6  the  specific  impulse  gain  realized  by  the  addition  of  PH3  to  the  fuel  is  graphed  as  a  function 
of  mass  per  cent  PH3  in  the  fuel  with  the  initial  gas  velocity  as  a  parameter.  These  calculations 
were  done  for  a  slightly  fuel-rich  reaction  mixture.  For  this  nozzle  geometry,  the  specific  impulse 
gain  appears  to  reach  a  maximum  at  2-3  per  cent  by  mass  phosphine  in  the  fuel.  While  the  values 
of  specific  impulse  gain  are  not  large,  no  effort  was  made  in  this  study  to  optimize  the  nozzle 
design  for  the  phosphine  fuel  additive. 
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Fig.  6  The  gain  in  specific  impulse  from  phosphine  addition.  The  phosphine  mass  per  cent  is 
indicated  on  the  abscissa. 
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